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ABSTRACT 

Absorption spectra of high redshift quasars suggest that the reionisation of cosmic hy- 
drogen was complete near z ~ 6. The dominant sources of ionising photons responsible 
for this reionisation are generally thought to be stars and quasars. In this paper we 
make a quantitative estimate of the relative contributions made by these sources. Our 
approach is to compute the evolution of the post overlap ionising background radia- 
tion by combining semi-analytic descriptions of reionisation in a clumpy medium with 
a model for the quasar luminosity function. Our overall model has two free parame- 
ters, the star formation efficiency and the minimum quasar luminosity. By adjusting 
these parameters, we constrain the relative con tributions made by stars and quasars 
through comparison with reported observations l)Fan et al.l 2005). We find that the rel- 
ative quasar contribution (at z = 5.7) to the ionising background was between 1.4% 
and 14.5%. The range of uncertainty is dominated by the unknown minimum quasar 
luminosity. 

Key words: cosmology: theory - galaxies: formation - intcrgalactic medium 



1 INTRODUCTION 

Observations of high redshift quasars imply that the reion- 
isation of cosmic hydrogen was completed by z ~ 6 
iFan et alJ fcOOSl hereafter F05); iGnedin fc Far] <2006ft : 
IWmTe^t^JTi| 2003)] . It is commonly believed that reioni- 
sation was d ue to UV photons prod uced by the first stars 
and quasars jBarkana fc Loebll20oJh . However, theoretical 
uncertainties and a lack of observational evidence make 
it difficult to estimate the relative contributions made by 
these sources. Although quasars and stars contribute equally 
to the ionising ba ckground by z ~ 3 IjKriss e t al. 2001; 
ISmette et al . 2002), it is often thought that stars dominated 
the UV flux at z > 6 (e.g. Fan et al. 2001). 

Attempts to make quantitative estimates of the contri- 
bution of quasars to the end of the reionisation era have 
proceeded o n several front s and led to inconsistent results. 
For example. iMadau et al] (119991 hereafter MHR99) estab- 
lish a threshold photon emission rate at z = 5 below which 
it would not be possible to maintain the reionisation of the 
inter-galactic medium (IGM). They then compute the pho- 
ton emission rate coinciding with the known QSO popula- 
tion at z ~ 5 and conclude that the quasar contribution to 
the budge t of io nising photons was ~ 10 -2 . Analogously, 
IFan et all ll200ll) compute the quasar contribution to the 
photon emission rate at z ~ 6 using the luminosity function 
(LF) computed from Sloan Digital Sky Survey (SDSS) ob- 
servations of high redshift quasars. They conclude that the 



quasar contribution to the ionising background at z — 6 was 
also ~ 10~ 2 and suggest that quasars could not be a major 
contributor of ionising photons at z ~ 6 without invoking a 
large negative luminos ity evolution in the LF. S imilar con- 
cl usions are drawn bvfYan fc Windhorst! ll2004l) and noted 
bv lStiavelli et all (120041) in support of the view that the ma- 
jor contribution to the ionising background at z ~ 6 comes 
fr om star forming galax ies. A different approach was taken 
bv lDiikstra et all (I2004T) , who show that the unresolved com- 
ponent in the observed soft X-ray background is beneath the 
level that would be contributed by a population of quasars 
that could fully ionise the universe at z ~ 6. 

On the other hand [Bunker et all ]2004l) show that the 
contribution to the ionising background from star forming 
galaxies identified in the Hubble Ultra-Deep Field (UDF) 
and Great Observatories Origins Deep Survey (GOODS) is 
insufficient to produce reionisation at z ~ 6. The missing 
ionising flux could be generated by quasars, or it could be 
generated by dwarf galaxies which w ere suppressed following 
reioni sation dWvithe fc Loebll2006l) . Alternatively. iMeiksinl 
(2005) finds that a quasar population described by a statis- 
tical fit to the double power law form of the LF dBovle et al] 
1988) falls short of the required flux to reionise the universe 
at z ~ 6 only by a factor of 2 — 3 and could therefore con- 
tribute this missing ionizing flux. 

This range of conclusions reflects the differing ap- 
proaches used to describe the ionising sources and more 
significantly, the different photo-emission rates assumed to 
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be necessary for reionisation. In the aforementioned analy- 
ses, reionisation in an inhomogeneous (clumpy) universe was 
treated using a clumping factor, C = (n )/ (n) (where n is 
the hydrogen number density), to account for the increased 
recombination rate. Larger values of C increase the emis- 
sion rate necessary to achieve reionisation. MHR99 consider 
a value of C = 30 correspo nding to the numerical simu- 
lation of iGnedin fc Qstrikerl jl997h . however they note the 
considerable unc e rtainty in the appl i cabilit y of this value. 
iFan et alJ fcOOll). lYan fc Windhorst] J2004h . IStiavelli et alJ 
ll2004) andlBunker et alJ l)2004ll use this same val ue for C. 
iDi ikstra et al.1 (2004) use a value of C = 10 and iMeiksinl 
J2005I) assumes values of 1 < C < 10. 

Much of the uncertainty in the ionising photon budget 
therefore centers around the value of the clumping factor. 
It is important to note that the clumping factor should be 
computed over the volume of IGM that must be maintained 
in an ionised state. This volume should not include the col- 
lapsed regions inside the virial radius of dark matter halo 
hosts containing star forming galaxies (the additional flux 
required to keep this overdense region ionised is accounted 
for by the value of escape fraction of ionising photons) . The 
volume over which the clumping factor should be computed 
also excludes those overdense regions in the IGM that form 
Lyman limit systems. These overdense clouds set the mean- 
free-path (MFP) for ionising photons and the recombina- 
tion rate should be computed only up to densities separat- 
ing the ionised IGM and these clouds. A framework within 
which the clumping factor appropriate for reionisation in 
an inhomogeneous IGM may be computed was described in 
iMiralda-Escude et alJ fcOOol hereafter M-E00). We follow 
the formalism developed in that work throughout this paper. 

We begin by describing our models for the reionisation 
history and quasar LF in § [5] We then review our proce- 
dures for calculation of the ionising background using these 
models in § |3] and present results from our analysis in § 01 
before discussing our conclusions in § |K| Throughout the 
paper we assume val ues for cosmological parameters based 
on WMAP3 results (ISpergel et al.ll2006ft but obtained by 
averaging over the observational data sets from WMAP3, 
SDSS, HST, SN Astier and BAO. The resulting parameters 
are fi A = .743, fi M = .257, fi B = .0437 and h = .718 (Renyue 
Cen, private communication). In computation of the mass 
function we assume a primordial power spectrum defined by 
a power law with index n = . 95, an exact transfer function 
given bv lBardeen et ail Jl986f) and rms mass density fluctu- 
ations with a sphere of radius Rg = 8h~ Mpc of as = 0.765. 



2 SEMI-ANALYTIC MODELS FOR 

REIONISATION AND THE QUASAR 
LUMINOSITY FUNCTION 

The density field in an inhomogeneous universe can be de- 
scribed by the overdensity, Aj = pi/~p~- M-E00 show how 
the effective recombination rate in an inhomogeneous uni- 
verse may be determined dynamically by consideration of 
the maximum overdensity (A cr it) penetrated by ionising 
photons within Hn regions (i.e. th e volume of the I GM to 
be maintained in an ionised state) . fWvithe fe Loebl ([2003a, 
hereafter WL03) employed this prescription for the recombi- 
nation rate into a semi-analytic model of reionisation. In this 



model Acrit is fixed prior to overlap, after which it evolves 
with redshift [see equation (6) and relevant discussion in 
WL03]. In this formalism it is possible to compute the ef- 
fective clumping factor in regions where Aj < A cr i t such 
that C — C[z, A cr it (2)]. For example, assuming overlap to 
occur at z ~ 6 and A cr it= 20 prior to overlap (following 
WL03), C(z = 6) ~ 2.3. Thus, the effective clumping fac- 
tor in ionised regions is considerably smaller than the value 
often assumed, which is over the whole IGM. 

The sources of ionising photons in this reionisation 
model were assumed to be quasars, in addition to popu- 
lation II (popll) and population III stars. In the current 
work we limit our attention to reionisation due to popll 
stars (Gnedin & Fan 2006), which govern the final stages of 
reionisation even in the presence of an earlier partial or full 
reionisation by population III stars (Wyithe & Loeb 2003b) . 
The quasar LF model of Wyithe & Loeb (2003b) success- 
fully reproduces 1 observations of the quasar number density 
over a large range of luminosities at z > 3. 

Our model includes two free parameters. First, the min- 
imum luminosity for quasars L m i u , which is not predicted 
by the model. Second, the parameter f^ BC represents the 
product of star formation efficiency and escape fraction, and 
characterises the contributions made by stars. Our approach 
is to compute a reionisation history given a particular set of 
L m i n and /* sc . With this history in place we then compute 
the evolution of the background radiation field due to these 
same sources. Post overlap, ionising photons will experience 
attenuation due to residual overdense pockets of Hi gas. We 
use the reionisation model of WL03 to describe the ionis- 
ing photon MFP, and subsequently derive the attenuation 
of ionising photons. We then compute the flux at the Lyman 
limit (y^ — 3.29 x 10 15 Hz) in the IGM due to sources imme- 
diate to each epoch, in addition to redshifted contributions 
from earlier epochs. 

We note that Hi Lyman limit photons (13.6 eV) are 
incapable of ionising helium (He) [24.6 eV (Hen), 54.4 eV 
(Hem)]. We therefore neglect He when computing the in- 
tensity of the ionising background. However, the presence of 
He is fully incorporated in the reionisation model of WL03 
which we use to describe the history of reionisation. 



3 EVALUATING FLUX AT THE LYMAN 
LIMIT 

In this section we review the process for calculating the ion- 
ising background flux at a particular redshift. We assume the 
background to be generated by a combination of stars and 
quasars. We first define zo to be the redshift at which the 
flux is to be evaluated. The flux at the Lyman limit is nor- 
malised in physical units of J21 (10 _21 ergs/sec/Hz/cm 2 /sr), 
and at redshift zq is related to the energy density by 



Jai(«o) 



1 



4tt dVdv 10- 21 



(1 + z f 



(1) 



1 This model is less successful at reproducing the number of very 
luminous quasars at z ~ 6 using the WMAP third year cosmology, 
indicating the necessity of revision of its physical basis. However 
we maintain the LF described in Wyithe & Loeb (2003b) as an 
empirically successful description. 
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where c is the speed of light and — dvlv °^ * s ^ e ener Sy P er 
unit frequency interval, per co-moving volume, at frequency 
v-l, and zo. Contributions to the radiation field at z/l (and 
zo) from sources at redshift z were emitted at frequency 
f z — }+ z z/ L . Note that v z remains below the He Lyman 
limit (24.6 eV) whilst < 2. The semi-analytic model of 
the reionisation history predicts an average redshift of over- 
lap, which we refer to as z \. While the observed redshift of 
overlap is subject to some scatter among different lines-of- 
sight, we make the approximation that the universe becomes 
transparent to Lya photons at z a \ along all lines-of-sight. 

The ionising background flux contains contributions 
from sources at zo in addition to redshifted flux from sources 
at higher redshift. To compute the ionising background 
flux we consider contributions from sources with redshifts 
zo < 2 < z oi, and therefore assume the complete attenu- 
ation of ionising photons emitted prior to the redshift of 
overlap. The total co-moving energy density at redshift zo 
is 



dVdu 



dz 



d E^(z) _ T ( ZjZQ ) 
dVdudt e 



dt 

dz ' 



(2) 



where the exponential term describes the attenuation of ion- 
ising photons (§ 13.31 between redshifts z and zo- In this 
equation one might expect to find factors of ^ j due to 
expansion. However, energy loss due to redshifting of the 
photon frequency is cancelled by the redshifting of the fre- 
quency interval. The term d d y V d 1^ represents the frequency 
dependent (co-moving) energy density per unit time gener- 
ated by ionising sources. Finally, the relation between proper 
time and redshift is 

dt 



— = [Ho(i + z )^/(i + zym M + nAr\ (3) 

where Ho is the current value of the Hubble parameter. 
Equations Q & J^J describe the ionising flux at zo given 
contributions of ionising sources at higher redshifts. In the 
next two sub-sections we describe the calculation of the 
emissivity of quasars and stars, which provide the sources 
of ionising radiation. 



3.1 Lyman limit photons from Quasars 

Our estimate of the emissivity of quasars is based on the B- 
band quasar LF [<6(Lb,z)] ; i.e. the number density per co- 
moving volume per solar B- band luminosity Lb- We use the 
theoretically derived LF of lWvithe fc Loebl 112003 a) , which 
successfully reproduces known properties of the observed 
quasar LF at both the bright and faint ends at z > 2.5, 
including the isolated observations discussed in § 14.31 This 
model is based on the physics of galaxy mergers, accretion 
powered luminosity and the self-regulated growth of super- 
massive black holes. The total B-band luminosity [Lb (•?)] 
per volume is calculated directly from &(Lb, z) 



dL B tot (z) 
dV 



dL B <E>(L B , z) x L B . 



(4) 



While this integral converges at high luminosities, the shal- 
low slope of the LF at low luminosities means that the value 
of the integral is sensitive to the lower limit L m i n , which is 



highly uncertain. We therefore treat L m i n as a free parame- 
ter in our description of quasar emissivity. To convert from 
Lb to the luminosity at the Lyman limit (L Vl ) we use the 
relation described bv lSchirber fc Bullock! (12003ft . 



7?„ L = — - = 10 ergs s Hz Le 
Lb 



(•») 



In order to compute the Lyman limit flux at zo, we require 
the contribution from intrinsically higher frequency photons 
from higher redshifts. We assume a powerlaw spectrum L„ oc 
v~ a , whence we obtain 

1 + 2 



^ = ^^ =r ^{i + ^) ■ (6) 

iTelfer et all <|2002T) find a = 1.57 in the wavelength range 
500 < A < 1200 A using (radio quiet) AGN spectra 
obtained from the H ST UV survey. We adopt this value 
which is also used by ISchirber fc Bullock! i2003l) in deter- 
mining 77„ L . In this work we compute the contribution of 
type-I quasars to reionisation. There may be additional ac- 
cretin g objects whose ionizing fl ux is blocked by a dusty 
torus IIElitzur fc Schlosmanl l2006'l. In unified schemes these 
obscured objects are the same as type-I quasars, however 
they are viewed from a different orientation. Our quasar 
LF accounts for the observed luminosity density of ion- 
ising radiation. We assume that all ionising photons es- 
cape the host galaxy [note however that the Lya halos 
around some QSOs indicate r eprocessing of this ionising flux 
jFrancis fc McDonnell l2006ty . Finally, using equation @, 
we find the energy density due to quasars at frequency z/ L 
and redshift zq, 



d 2 El^{zo) 
dVdv 



?7" L / dz 



° dL B tot (z) (l + zo 



dV 



l+Z 



dt 

dz 



(7) 



This energy density may be converted to a flux and ex- 
pressed in units of J21 using equation 



3.2 Lyman limit photons from Stars 

We next describe the contribution made by popll stars. To 
begin, we describe the spectral energy distribution (SED) 
of popll star forming galaxies using the model presented in 
iLeitherer etafl lll999l) . The SED has units of ff"-, , where 

1 — ' * " ■ dvdtdM ' 

M is expressed in units of (baryonic) solar masses per year. 
In the instance that ionising photons are produced primarily 
in starbursts, with lifetimes much shorter than the Hubble 
time, we may express the star formation rate per unit time 

as 



dM , . 

■dv {z) 



r dF(z) 
f dt Ph ' 

111 vear 



(8) 



where pb is the co-moving baryonic mass density, and F(z) 
is the collapsed fraction of mass in halos above a critical 
mass at z. The factor /* (star formation efficiency) de- 
scribes the fraction of collapsed matter that participates in 
star formation. This f raction is largely unknown, however 
IWvithe fc Loebl J2006I) constrain /* to be ~ 10 - 15% by 
fitting to the galaxy LF at high redshift. 
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To compute F(z) we use the Press- Schechter (1976) 
mass function with the modification of ISheth fc Tormenl 
( 1999) to describe the evolution of the collapsed fraction 
[F(z)] within halos above a critical virial temperature. In 
a cold neutral IGM beyond the redshift of reionisation, 
this critical virial temperature is set by the temperature 
(Tn ~ 10 4 K) above which efficient atomic hydrogen cool- 
ing promotes star formation. Following the reionisation of a 
region, the Jeans mass in the he ated IGM limits accretio n 
to halos above Ti ~ 2.5 x 10 5 K jThoul fc Weinberdll99ri . 
The correspondi ng virialised mass can b e determined using 
equation (26) in Barka na fc Loebl i200lTl . 

We may therefore write the time derivative of the col- 
lapsed fraction 



dF 



<lty 



dz 



dF(z,T!) 
dz 



+ [1 - Q m (z)] 



dF(z,T N ) 



dz 



dty 



(9) 



here Q m is the ionised baryonic mass fraction in the universe. 

To describe the ionising flux from stars we require one 
further parameter. We have assumed that all ionising pho- 
tons produced by quasars escape the host galaxy (§ 13.11 . 
However, due to the softer spectra of galaxies it is likely 
that only a fraction of ionising photons produced by stars 
enter the IGM. Therefore, an additional factor of / esc (the 
escape fraction) m ust be included when com puting the emis- 
sivity due to stars. Iclardi fc Ferraral i2005l) include a review 
of existing constraints on / csc which suggests that its value 
is < 15%. The star formation efficiency and escape fraction 
may be combined into a single free parameter (/e SC ) to de- 
scribe the contribution of stars in our model. 

Finally the energy density due to popll stars at zo may 

be computed using equation [5] with d 

d 3 Ev z {z) _ d z E Vz {z) dF(z) 



■X") 



dVdvdt 



dvdtdM dt 



Ph /esc- (10) 



This energy density may then be converted to a flux and 
combined with the contribution due to quasars. 



the MFP is set by residual Hi in ionisation equilibrium with 
the background flux as well as Lyman limit systems. Our 
estimate of the MFP therefore provides an upper limit. Fur- 
thermore, our reionisation model assumes ionising photons 
to be instantly absorbed, an assumption which becomes less 
applicable post overlap. 

Finally the MFP as a function of redshift may be used 
to compute the attenuation of ionising photons between red- 
shifts z and zq. The resulting optical depth is 



t(z, Zo) = 



' cdt 1 

dz-j—T-- 

dz m 



(13) 



4 RESULTS 

The model described in the previous section has two free 
parameters L m i n and /^ c . For any combination of these pa- 
rameters we are able to compute a reionisation history, the 
evolution of the MFP, the evolution of the ionising back- 
ground, and the relative contributions of stars and quasars 
to the ionising background. In this section we first sum- 
marise the observations to which we compare our model be- 
fore describing the results of this comparison. 



4.1 Observational estimates of the ionising 
background and mean-free-path 

F05 have analyzed the absorption spectra of 19 high redshift 
quasars. Based on this analysis F05 present estimates of the 
neutral fraction at a range of redshifts near the end of the 
reionisation era, as well as the background ionising flux and 
MFP. In this paper it is our aim to constrain the contribu- 
tions of quasars and stars to the ionising background based 
on these observations. 

The observed estimates of the ionising background are 
quoted as photoionisation rates (F05), derived from the ob- 
served (Lya) optical depth. These photoionisation rates may 
be related to ionising background fluxes using 



3.3 Attenuation of ionising photons 

Given the luminosity density we may find the ionising back- 
ground using equation (2), following calculation of the op- 
tical depth (t) and thereby the attenuation of ionising pho- 
tons. Following overlap, we use the approach of M-E00 to 
estimate the MFP (Aj) of ionising photons as a function of 
redshift (zt). 

Xi — Ao(l — F v )~ 2 ^ 3 . (11) 

Here we use the reionisation model of WL03 to compute the 
redshift evolution of the the volume filling factor in ionised 
regions. 

/•A cr it(z) 

F v = / dAP(A), (12) 

Jo 

where [P(A)] is the volume weighted probability distri- 
bution of the overdensity (M-E00). The product XoH = 
60kms~ 1 was obtained from comparison to the scales of Lya 
forest structures in simulations at z = 3 (M-E00). In reality, 



r = 4tt / -?-<T v dv, (14) 

where cr„ is the Hi photoionisation cross-section and h p is 
Planck's constant. Using equation 1141 and the spectra de- 
scribing our sources, we find Ti2 = I.69J21 (where Fi2 is the 
ionising rate in units of 10 _12 s _1 ) for a background in which 
stars are the dominant source. Similarly, if the background 
radiation is dominated by quasars we obtain Fi2 = 2.64J21. 
In comparing our model to observation we normalise the ob- 
served flux using a weighted sum of the instantaneous flux 
contributions from each of the stellar and quasar sources. 

F05 also estimate the ionising photon MFP. They com- 
pute MFPs using the frequency averaged cross-section ((ov)) 
to Lyman limit absorption assuming an ionising background 
spectrum of the form J„ oc v~ a . They compute (o v ) in the 
frequency range between the Hi and the Hei Lyman limits, 
and assume a = 5 which describes an ion ising background 
dominated by stars (Bark ana fc LoetfeOOlT) . In § 13.21 we out- 
li ned our calculation of s tellar emissivity using the spectrum 
of lLeitherer et aD fl999) . We therefore use this spectrum to 
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redshift redshift 

Figure 1. The best fit model (see text) corresponding to f* BC = .00815, L m in= 10 12 ' 5 Lb and A cr it=20. (a) The reionisation history. 
Solid, dotted and dashed lines represent the evolution of the volume and mass filling fractions [F v (Hn),F m (Hn),.F v (Heni)]. (b) The 
fractional contribution of quasars to the instantaneous production rate of ionising photons (solid line). The fractional contribution of 
quasars to the cumulative number of ionising photons produced by z (dashed line), (c) The evolution of the MFP. Our model (solid line) 
and observational data points from F05. (d) The redshift evolution of J21. The model combined flux (solid line) and contributions from 
stellar sources (dotted line) and quasars (dashed line). The data points are from F05. 



adjust the values of (o v ) from F05 (assuming stellar domi- 
nance) and present data points using this modified value. 

4.2 Comparison with observation 

In this section we show results for a particular model in 
which we assume L m i n = 10 12 ' 5 Lb - In this case we find that 
the model best fits the ionising background observations for 
the value of /o SC = -00815. For these parameters we find that 
model fit to observations of the ionising background yields a 
reduced x 2 value of 1.3. In a subsequent section, we explore 
the range of values for L m i n permitted by observation. 

For our representative model, Figure 1 shows the red- 
shift evolution of four quantities computed. Here we have 
assumed a critical overdensity of A cr it = 20 prior to overlap, 
which occurs in this case at z a \ ~ 6. Figure 1(a) shows the 
full reionisation history. The volume fraction of Hn (solid 
line), the mass fraction of Hn (dotted line) and the vol- 
ume filling fraction of doubly ionised Hem (dashed line) 
are shown. Note that overlap is marked by the point when 
the volume filling fraction of Hn reaches unity. Figure 1(b) 
shows the resulting fractional contribution of quasars to the 
instantaneous production rate (solid line) of ionising pho- 



tons. The dashed line shows the fractional contribution of 
quasars to the cumulative number of ionising photons pro- 
duced by redshift z. Ionising photons above the Hei Lyman 
limit also contribute to the reionisation of Hi (WL03). We 
therefore consider the contribution of ionising photons made 
by quasars and stars at all frequencies above the Hi Lyman 
limit. The contribution of quasars to the instantaneous pro- 
duction rate rises towards low redshift and becomes ~ 40% 
that of stars by z ~ 2.5, in broad agreeme nt with obser- 
vation llKriss et alJl200ltlSmette et alj l2002). Cumulatively 
the stellar contribution remains dominant even at low red- 
shifts. 

In Figure 1(c) we show the evolution of the MFP (solid 
line) with comparison to the points from F05. It is clear from 
this plot that the models abilty to predict the MFP dimin- 
ishes towards lower redshift. In Figure 1(d), we plot J21 as a 
function of redshift. The evolution of the total (model) flux 
due to quasars and stars is shown as a solid line. The con- 
tribution to the ionising background due to quasars (stars) 
alone is shown as a dashed line (dotted line). In Figures (c) 
and (d), excepting the two highest redshift bins the data 
points show the mean value for J21 with 1 — a error bars. 
In the two highest redshift bins optical depth measurements 
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Figure 2. The a posteriori probability (P) assuming values for X/ m i n (X/B) = 10 11 (solid), 10 1175 (dashed), 10 12,5 (dot-dashed), (a) 



The cumulative distribution as a function of /* 
production rate at z = 5.7. 



(b) dP/dF where F is the fractional quasar contribution to the instantaneous photon 



were not always possible due to complete Gunn-Peterson 
troughs along some lines of sight (see F05). In these cases 
the mean and standard deviation were calculated using lower 
limits on the optical depth. Hence these data points repre- 
sent upper limits on J21. 

4.3 Constraining f^ sc and the quasar flux 
contribution to the IGM 

In the previous sub-section we discussed the results from a 
particular model which assumed a value for L m i n = 10 12 ' 5 Lb. 
In this sub-section we investigate values of /* sc , given various 
minimum quasar luminosities. In particular, for a range of 
Lmin we compute the a posteriori probability of our model 
flux, as a function of the free parameter /* sc , given the like- 
lihoods of our model fits. We include an a priori probability 
distribution which is flat in the logarithm and consider a 
range for / c * sc 6 [.006, .01]. 

At low redshift (z<.5) quasars are observed at lumi- 
nosities corre sponding to ~ 10 7 L B iBovle et alj EoOO). 
iBarger et alJ (120031) have inferred from X-ray observations 
a quasar with luminosity of ~ 10 113 Lb at z = 5.7. At 
z = 5.85 observations have revealed the presence of quasars 
at lu minosities corresponding to ~ 10 12 ' 8 Lb dCool et alJ 
2006). This suggests an upper limit to L m i n as even at 
high redshift we can see t hese quasars directly. At z ~ 4.5 
iDiikstra fc Wvithei ll2006l) find a downturn in the density 
of AGN with L < 10 n L B , based on the absence of AGN in 
Lya surveys. However, current data is insufficient to strongly 
constrain the lower limit of QSO activity. Therefore, within 
this range we construct conditional probability distributions 
at Lmin (Lb) = 10 11 , 10 1175 , 10 125 . In each case probability 
distributions are extracted from the aforementioned likeli- 
hood analysis. For display we normalise the probability den- 
sity functions such that the maximum likelihood is set to 
one. 



In Figure[5Ja) we show the resulting cumulative a poste- 
riori conditional probability distribution for /e SC . The mod- 
els which best fit the ionising background observations, for 
each given L m t n , have corresponding values for /* sc in the 
narrow interval .0074 < /e SC < .008. Furthermore, for each 
particular L m i n the distribution is restricted to a narrow in- 
terval in /* sc space. This is a reflection of the sensivity of the 
goodness of fit to /o SC . The best fits for each L m i n correspond 
to overlap redshifts suggested by the reionisation model at 
z ~ 6. This is slightly lower than the recent numerical anal- 
ysis of Gnedin & Fan (2006) which suggested overlap occurs 
between 6.1 < z a \ < 6.3. Better fits are obtained for larger 
values of L m i n . However, this reflects systematic uncertainty 
in the model rather than a real constraint. In Figure[5f b) we 
show the corresponding distribution 4£. — -jp— jjp , where 
F represents the fractional contribution of quasars to the 
instantaneous photon production rate at z = 5.7. For the 
values of L m i n assumed in Figure |5|the most probable (frac- 
tional) quasar contributions to the ionising background are 
F ~ .145, .061, .014. This range covers the systematic uncer- 
tainty introduced by L m i n . 

The models presented thus far have assumed a value of 
A cr it = 20 beyond z i. We have repeated our analysis for 
different choices of A cr i t and find that while the results are 
similar, the fit to J21 improves as A cr i t is increased between 
5 and 20. In particular, for higher A cr i t the flux immedi- 
ately following overlap rises more sharply due to both the 
increased /* sc that is required to match observation and also 
the higher value of Fm at z \. 

4.4 Additional Uncertainties 

Our model contains several additional uncertainties. These 
include, but are not limited to, the following. First, the 
model LF describes the density of high redshift quasars. The 
uncertainty in the overall level of quasar flux at z ~ 6 (where 
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there is only 1 point) is at least a factor of 2. This uncer- 
tainty needs to be added to the systematic uncertainty in 
Lmin. Second, the MFP computed by our model is an upper 
limit. Indeed comparison of our model with the results of 
Fan et al. (2005) suggest an overestimate by a factor of ~ 2 
at z < 5.5. This overestimate results in an underestimate 
of the attenuation in our model. We may therefore under- 
estimate the quantity /* sc required for a good description 
of the ionising background, and therefore overestimate the 
fractional contribution to quasars. Thirdly, our model does 
not consider the possible evolution of L m j n with redshift. 



5 CONCLUSIONS 

In this paper we have estimated the relative contributions 
made by quasars and stars to the ionising background at 
the end of the reionisation epoch. We have characterised 
these contributions by two free parameters. L m m which de- 
fines the minimum luminosity of contributing quasars and 
/esc which describes the product of star formation efficiency 
and escape fraction of ionising photons from galaxies. We 
compute the relative quasar contribution for the range of 
values, L m in(iB) = 10 11 , 10 11 ' 75 , 10 12 ' 5 . We find the relative 
contributions made by quasars to the ionising background 
at z = 5.7 are F ~ .145, .061, .014, respectively. The value of 
Lmhi provides the greatest uncertainty. 

We are able to quantitatively constrain the contribu- 
tion of quasars to the ionising background because our model 
(following M-E00) computes the clumping factor in the IGM 
which has represented a major uncertainty in previous stud- 
ies. There are a few independent checks which enhance our 
confidence in the results. Firstly, our quasar LF accurately 
reproduces observations at z > 2.5. This quasar popula- 
tion was responsible for the double reionisation of helium, 
which for values of Z/ m i n (-I/B) ~ 10 11 is predicted by our 
reionisation mod el to occur at 3 < z < 3.5, in agreemen t 
with observation feernardi et alJl2003t TTheuns et alJl2002l) . 
Secondly, the contribution of quasars to the ionising back- 
ground becomes c omparable to that of stars by z ~ 2.5 
(Kriss et al]l200ll : ISmette et al.ll2002h . Thirdly, our model 
predicts /e SC = .00815 which is consistent with the pr oduct 
of the values f* ~ 10 - 15% jWvithe fc LoedEiolj) and 
fesc < 15%. ijCiardi fc Ferraral 12005ft which are estimated 
from other studies. In future, improvements to the analy- 
sis presented here could be made by adding a population 
of quasars to cosmological simulations of reionisation (e.g. 
Gnedin & Fan 2006) which more accurately reproduce the 
evolution of the ionising background. 
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